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The quantum walk (QW) is the term given to a family of algorithms governing the evolution of 
a discrete quantum system and as such has a founding role in the study of quantum computation. 
We contribute to the investigation of QW phenomena by performing a detailed numerical study of 
discrete-time quantum walks. In one dimension (ID), we compute the structure of the probability 
distribution, which is not a smooth curve but shows oscillatory features on all length scales. By 
analyzing walks up to = 1000000 steps, we discuss the scaling characteristics and limiting forms 
of the QW in both real and Fourier space. In 2D, with a view to ready experimental realization, we 
consider two types of QW, one based on a four-faced coin and the other on sequential flipping of 
a single two-faced coin. Both QWs may be generated using two two-faced coins, which in the first 
case are completely unentangled and in the second are maximally entangled. We draw on our ID 
results to characterize the properties of both walks, demonstrating maximal speed-up and emerging 
semi-classical behavior in the maximally entangled QW. 

PACS numbers: 03.67.Ac, 05.40.Fb 


I. INTRODUCTION 

Quantum walk(QW) is an extension of the classical 
random walk, in which the walker is controlled by a quan¬ 
tum mechanical variable [l|, 0 . The quantum interference 
from different walk paths brings about QW novel and in¬ 
teresting features. An motivation for work on QW is 
developing new quantum algorithms aiming at an 

exponential speed-up compared to classical algorithms 
for a certain classes of problem 

On the physical side, QW is a valuable model for 
studying evolution process from simple quantum proto¬ 
col: discrete-time QW @ and continuous-time QW in 
one dimension (ID) M-' Deep studies involve multiple 
coins m , multiple walkers [l^ , multiple dimensions [l^ , 
and multiple coins with different degrees of entanglement 
[11®. A recent resurgence of interest in the QW has 
come about to a significant extent because experimental 
capability has caught up with many of these theoreti¬ 
cal proposals, particularly to create walks with complex 
coin protocols and in two phase-space dimensions |17l - l^ . 
However, all of these QWs implementation are restricted 
to rather small numbers of steps, and some have inherent 
limits to the size of their phase space. 

Despite the large volume of work performed on the 
QW, we have not been able to find any that address the 
properties of the system at long evolution times, meaning 
at large values of the step number N. As the challenges 
in quantum computing move toward the large scale, it 
would appear that there is a need to understand the be¬ 
havior of quantum algorithms in the large-regime. The 
aim of this paper is to analyze QWs at large N to es¬ 
tablish their properties, both universal and specific, and 
thus effectively to gauge the flow and concentration of 
information in one particular set of algorithms. We will 


consider in detail the ID QW, where we have performed 
calculations up to A = 1000000 steps, to establish the 
scaling characteristics, spatial information content, uni¬ 
versality, and limiting form of the probability distribu¬ 
tion. We will then turn to different types of possible 
QW in 2D, where we use our ID knowledge both directly 
and to benchmark the additional forms of behavior that 
emerge, particularly entanglement, accelerated diffusion, 
and semi-classical limiting distributions. 

In this paper we perform a detailed analysis of the 
discrete-time ID QW and two fundamental 2D QWs. By 
calculating probabilities for large numbers A of steps, 
we investigate the destructive interference, the scaling 
properties, the frequency content, and the combination 
of QWs. The structure of the article is as follows. In 
Sec. II we review the classical random walk and a simple, 
symmetrical model for understanding the discrete-time 
QW in ID, including its analytical solution. Section III 
presents our numerical results for ID QWs up to large 
A and their analysis in both real and Fourier space. In 
Sec. IV we discuss two typical 2D QWs as a theoreti¬ 
cal basis of Sec. V. Borrowing from the understanding 
developed in ID, in Sec. V we provide the complete nu¬ 
merical characterization of these two 2D QWs. Section 
VI provides a brief summary. 


II. ONE-DIMENSIONAL RANDOM WALKS 

A. Classical Random Walk 

A walker standing at the origin of a line flips an un¬ 
biased coin and steps to the right if a head comes up or 
to the left if the result is a tail. After many flips, and 
taking a fixed number of steps based the coin state, the 
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FIG. 1. Probability distribution for a ID classical random 
walk with N = 100. 



FIG. 2. Probability distribution for a ID quantum walk with 
N = 100. 


location of the walker is an unknown, random position, 
but the probability distribution of this position is a def¬ 
inite statistical quantity. After N steps, the probability 
of the walker being at position x is 


the quantum system is described by the wave function 

N 

Hn) = (oil i) + h\ t))|i), (2) 
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i.e. this process, the discrete-time ID classical random 
walk, follows a binomial distribution. The walker’s posi¬ 
tion X may only be an even (odd) integer when JV is even 
(odd). In the continuum limit, which is approached for 
sufficiently large N, the distribution is Gaussian (Fig.[T|). 
Its standard deviation, which represents a mean propa¬ 
gation distance in a sample of many walkers, is tr = \/]V. 


B. Quantum Walks 


A QW refers to a random walk effected using a quan¬ 
tum device, and as such requiring a quantum mechan¬ 
ical description.Following the classical formulation of a 
coin and a walker, the QW differs fundamentally from 
its classical counterpart in that coin and walker are “en¬ 
tangled” in the same quantum “particle”. The essential 
consequence of this entanglement is that the propagation 
of the quantum walker depends not on its probability 
(Sec. IIA) but on its amplitude. This amplitude is sub¬ 
ject to quantum mechanical interference, which can be 
constructive or destructive, leading to completely uncon¬ 
ventional forms of probability distribution (Fig. [2]). 

The numerical results in Sec lIIIl based on following cal¬ 
culation. While there are many ways of describing such 
a walk (Ref. [ij for a review), here we introduce only the 
case of a discrete time symmetric QW on an open line. 
With a view to later application, we denote the two in¬ 
ternal “coin” states of the walker, or particle, as f and j,. 
In a classical walk the coin states are completely separate 
(t or 4 , with probabilities 0 or 1), whereas a quantum coin 
can occupy any superposition state a| f) + ^| i)- Thus 


where N is number of steps in the walk, i is the position 
index, and |i) the corresponding state. 

The QW is described by the evolution of this wave 
function under the quantum operation for successive 
steps. The probability distribution for finding the walker 
at position i (state |j)) is found from the trace over the 
coin states to be Pi = of + b‘f, with Pj = 1. 

The evolution has two substeps: (1) A rotation in the 
coin space, represented by a unitary operator U, which 
is often taken as the Hadamard transformation 


Uh 
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(2) A coin-dependent translation of the walker, described 
by the shift operator 


5 = I t)(t I ® E I*+i)(*i +1 i)a I ® E I* - (4) 

i i 


The coefficients ag and bg of the initial coin state are 
arbitrary, which has a profound influence on the sym¬ 
metry of the probability distribution of the QW [I|. To 
obtain a distribution symmetric under the application of 
the Hadamard operator, the initial coin state is written 
as 


IV’o) = ;^(lt)+iU))G|0), (5) 

A full step of the QW is 

IV^i) = su^l^g) = ^1 t)|l) + ^1 i)l - 1) (6) 

The probability distribution for a QW of A^ = 100 
steps is shown in Fig. [51 It is manifestly completely dif¬ 
ferent from the classical random walk (Fig.[T]), with clear 
maxima at values of i close to ± 0.7N (revealed at larger 
N to be iLN/\/2). The probability for I = 0 is close 
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to 0 , indicating that the origin of this counter-intuitive 
behavior is an almost complete destructive interference 
among the paths of the quantum walker returning to the 
origin. The standard deviation of the probability distri¬ 
bution in this QW is cr oc iV, indicating a linear spreading 
rate, which is one of the most important attributes of this 
evolution algorithm [Eq. ([5])]. We present a quantitative 
analysis of the properties of the QRW in Sec. III. The 
QW evolution process acts to propagate the quantum 
mechanical amplitudes, preserving the complete informa¬ 
tion content of the internal states. 

There are another two kinds of analytical descriptions 
do benefit to understand the evolution of QW. The differ¬ 
ence between the classical and quantum walks can be un¬ 
derstood from the non-commuting nature of the (matrix) 
quantum operators and its consequences for the interfer¬ 
ence of different walker paths. The Hadamard operator 
can be decomposed [13 as Uh = P + Q with 




and Q 




It is easy to see that P determines motion of the walker 
to the right and Q to the left. Evolution under the QW 
for N steps is represented as = {P + Q)^■ In the 
classical random walk one has = {p + q)^, where 
p = q = ^ are real numbers representing probabilities. 

One may conclude that the non-commutativity of the 
quantum operators determining different paths to the 
same walker position, encodes the interference of am¬ 
plitudes leading to the entirely unconventional quantum 
phenomena reflected in Figs. [T] and [2] We expand this 
method to study a novel 2D-QW behaviour in Sec IVBI 

We pay attention to another analytical solution^ in 
order to introduce the concept of Fourier analysis. The 
evolution operator has a more concise form in fc-space 
(Fourier space) than in x-space (real space), so the initial 
wave function may be transformed to and evolved in k- 
space. By inverse Fourier transformation, the real-space 
wave functions for the | f) and | ),) internal states are 
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and finally the probability of finding the walker at a given 
position X after a walk of t steps is given by 


P{x,t) = |^/>t(a:,t)p -f |^/> 4 ,(a;,t)|^ 


( 8 ) 


Eqs. (O and (|5]) contain no dynamics, because the term 
provides no more than a compact notation for the 
combination of N with fc, the spatial Fourier variable 
conjugate to the actual walker displacement x. 

A numerical calculation of the exact analytical solu¬ 
tion, contained in Eqs. © and ®, is shown as the red 
curves and points in Fig. [31 where we compute the prob¬ 
ability distribution in both real space [Fig. |31[a)] and 
Fourier space [Fig.jSKb)] for N = 100 and 1000. The blue 
curves show our numerical calculations based on method 
in Eqs. ([HI) with the initial state specified in Eq. ([5]). The 
results are identical up to a relative error of 10 “^ caused 
by the numerical integration of Eq. 0 . We comment in 
detail on the forms of these distributions in Sec. III. 

The importance of the values ±N/^/2, noted in 
Sec. IIB, is clearly evident in Fig. |3Ka), and it was de¬ 
duced in Ref. [l^l that the limiting distribution is con¬ 


centrated in the interval 


as —>■ oo. We will 


N N 

V2\ 

qualify this statement in Sec. IIlA. Although the ana¬ 
lytical solution @,113 gives the exact probabilities for 
any position and number of steps, in fact the expressions 
in Eq. 0 , which require numerical integration over com¬ 
plex quantities, are not easy to compute when N becomes 


large. In this regime, direct numerical calculation of the 
probability distribution is more straightforward, and we 
use this approach in Sec. Ill to reveal the properties and 
structure of the ID QW at large N. 


III. NUMERICAL STUDIES OF THE ID QW 


We now consider in detail the properties of the ID QW 
introduced in Sec. II. The probability distribution ob¬ 
tained from this quantum algorithm has a number of very 
striking properties in both real and Fourier space (Fig. [3]). 
Among them are the clear importance of ± N/\f2 noted 
above, the twin-peaked “envelope function” of the P{x) 
distributions with its remarkable zone of destructive in¬ 
terference around a; = 0 , the width of these peaks, and 
the rapid oscillation of the functions (both P{x) and 
F{k)) at high spatial frequencies within the envelope. 
We clarify immediately that this oscillatory behavior is 
“real” in the sense that Figs. |3| and |3| show only the prob¬ 
abilities at even steps 0, ±2, ±4, ..., with the zero- 
probability odd steps not shown; the QRW contains ad¬ 
ditional oscillations in space between the “microscale” of 
individual steps and the “macroscale” of the walk length. 
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FIG. 3. (color online) Probability distribution for the ID QW, 
comparing the exact expression of Eqs. 0 and 0, shown 
in red, with numerical calculations based on Eqs. © and 
© shown in blue, for walks of = 100 and 1000 steps, 
(a) Probability distribution P{x) in real space, (b) Fourier 
components F{k) of P{x). Insets for N = 1000 show (a) the 
probability oscillations near a; = 0 and (b) beats in the Fourier 
envelope near fc = tt (see text). 



FIG. 4. Numerical results for the ID QW with different num¬ 
bers N of steps, scaled to the same horizontal axis. Proba¬ 
bilities are shown only for even values of the site position x, 
with zero values on odd sites excluded from the curves. 
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A. Nl\pl Property 

We begin by analyzing the most obvious feature of the 
QW, which is the tendency for the probability distribu¬ 
tion to peak around 0.7N. From Sec. IIC it is clear that 
the factor 1 / plays an important role in the analytic 
solution and in physical terms it would appear to mark 
the crossover in behavior from a low probability arising 
due to destructive interference to a low probability aris¬ 
ing simply from the extremely low likelihood of having 
more than 85% of the steps of an unbiased random walk 
be in the same direction. In Fig. 2] we show numeri¬ 
cal results for the probability distributions of ID QWs 
with four different values of N, with the position axis 
normalized by A. As N increases, the distributions ex¬ 
hibit both increasingly oscillatory behavior, which we an¬ 
alyze in Sec. IIIB, and a peaking of the envelope function, 
which becomes sharper as it converges towards a maxi- 


FIG. 5. (color online) Scaling of the position Xmax for maxi¬ 
mum probability with total step number N in the ID QW. 

mum probability close to step N(\/2 (Sec. IIIC). 

This convergence in shown in more detail for values 
up to iV = 1000000 in Fig. © The macroscopic feature 
is indeed a convergence towards A/\/2. Further, on a 
relative scale the distribution appears to tend towards 15- 
functions centered at ±A/-\/2. However, we caution that 
this is not the complete story and we consider the asym¬ 
metric envelope shape in Sec. IIIC. The point of maxi¬ 
mum probability is in fact Xmax = 70684 for N = 100000 
and Xjnax = 707050 for N = 1000000, while the exact 
value of l/-\/2 is 0.707106. Thus in fact there are still 
more than 25 even steps of finite probability separating 
a^max from iV/-\/2 for N = 1000000 [see Figs. [7l[i) and 
©j)]. The probability at a; = N/y/2, shown in Fig.© is 
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FIG. 6. Representative probability data from different points 
on the distribntion, x = 0, x = Af/2%/2, x = N/2, x = Xmax, 
X = N/y/2, and x = 0.7072A'^, for ID QWs of step lengths 
from N — 1000 to = 1000000. Where relevant, all data are 
taken from the upper envelope of the distribution. 


always close to one half of P(xma.x), and this point marks 
the approximate crossover where P{x,N) changes from 
algebraic to exponential decay with N (below). There are 
always points of finite P{x) beyond x = 
more exact statement of the result of Ref. is that the 

normalized support converges to the interval ~ 

as iV —>■ oo. The probability of a walker passing beyond 
this interval vanishes more rapidly than interference ef¬ 
fects can cause it to grow, and the net consequence of the 
destructive interference between paths is to “pile up” the 
probability close to (but mostly below) the limits of the 
interval. 

Further insight into the nature of the (upper) envelope 
function may be obtained by considering the probabili¬ 
ties at different representative positions on the normal¬ 
ized x-axis, as shown in Fig. [51 The probabilities Pn{0), 
P/v(7V/2-\/2), and Pn{N/2) all fall linearly with 1/fV, sug¬ 
gesting a constant weight if binned into intervals whose 
width scales with N. It is worth noting that the weight at 
position zero, which has the maximum number of inter¬ 
fering paths, does not vanish completely due to destruc¬ 
tive interference in any finite-length QW. By contrast, 
the probability at positions x = x^ax and x = N/\/2 
scales as P oc A^“2/3 [specifically, P{xniax) = 
and P{N/\/2) = 0.44P(a:niax)], accounting for the sharp¬ 
ening of the distribution peaks with increasing N. We 
return to the question of the peak shape in Sec. IIIC. 
As noted in the preceding paragraph, beyond x = N/\/2 
there is a very abrupt change in the form of P{N) to 
an exponential decay, as shown in Fig. for the point 
X = 0.7072A^ and in the sudden loss of oscillatory behav¬ 
ior in Figs.[7Ki) and[7Kj); this we will also analyze in more 
detail in Sec. IIIC. 


B. Oscillatory Behavior 

We turn next to the question of the oscillatory behavior 
of the probability distribution within its envelope func¬ 
tion. We stress again that this has nothing to do with 
the period-2 oscillation created by the fact that walkers 
alternate between odd and even sites at successive steps 
of the walk. We begin by showing in Fig. 0 the qualita¬ 
tive nature of the oscillations in the real-space probability 
distribution function for selected regions of the interval, 
using different values of N to highlight their universal 
nature. 

One of the most important properties of the oscil¬ 
lations is that their effective wavelength, in terms of 
the fundamental step size, appears to decrease towards 
larger values of a;. At the center of the distribution, 
as shown in Figs. (TJa) and EKb), they have relatively 
long wavelengths, but towards the edges these become 
shorter [Figs.[7Kc) andlT^d)] until in the region x ~ 0.6A^ 
they have only twice the fundamental length [Figs. [T^g) 
and[7|[h)]. Towards the center of each half of the dis¬ 
tribution, some beat-like structures develop [Figs. [71[e) 
and[7l[f)]. As the peak is approached, the frequencies 
drop and the oscillations vanish suddenly at x = N/ 
[Sec. IIIA, Figs.[7l[i) and[7l[j)]. We draw attention to the 
fact that the distribution also has an effective lower en¬ 
velope function, in that the probability is never zero on 
any even points and in fact is very much larger than the 
size of the oscillations close to x = 0 [Figs.[7l[a) and[7l[b)], 
but we do not analyze this further here. 

To quantify the nature of the oscillations, we count the 
numbers of peaks in the probability distribution within 
different intervals and for QWs of different N. First of all, 
by counting the total number of peaks it is clear (Fig. |S]) 
that this is linearly proportional to (approximately 8.5% 
of) N. Subdividing the QW into intervals of fixed length 
and counting the peaks in each of these gives the re¬ 
sults shown in Table HI By comparing these peak counts 
horizontally, meaning for different intervals within walks 
of the same N, a steady increase in frequency becomes 
apparent out to x ~ 0.6A^, from very long-wavelength 
(periods of 50 or more steps) oscillations near x = 0 
to extremely rapid (period-4) ones when x is a signif- 
cant fraction of l/-\/2- We remind the reader here that 
a period of 4 in a system where odd sites have probabil¬ 
ity zero is essentially a max-min-max-min-... structure 
within the envelope [Figs.[71[g) and[71[h)]. Thus spatial in¬ 
formation about the QW is truly contained on all length 
scales. By comparing the peak counts vertically, mean¬ 
ing for different values of A^, it becomes apparent that 
corresponding regions have precisely the same frequen¬ 
cies, with the maximum frequency occurring in the re¬ 
gion around x = 0.58N. Thus the spatial modulation of 
the QW is a quantity independent of N; although QWs 
of different N cannot be called self-similar, they do share 
similarities in particular aspects. 

Other forms of similarity and scaling appear in particu¬ 
lar segments of the probability distribution. Considering 
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TABLE I. Numbers of peaks in the probability distribution P{x) found within windows of width 100 steps for QWs of = 1000, 
10000, and 100000. 


interval {N = 1000) 

[0,100] 

[100,200] 

[200,300] 

[300,400] 

[400,500] 

[500,600] 

[600,700] 

number of peaks 

2 

5 

8 

12 

17 

23 

17 

interval {N = 10000) 

[500,600] 

[1500,1600] 

[2500,2600] 

[3500,3600] 

[4500,4600] 

[5500,5600] 

[6500,6600] 

number of peaks 

2 

5 

8 

12 

17 

23 

17 

interval {N — 100000) 

[500,600] 

[15500,15600] 

[25500,25600] 

[35500,35600] 

[45500,45600] 

[55500,55600] 

[65500,65600] 

number of peaks 

2 

5 

8 

12 

17 

23 

16 

interval {N = 10000) 

[5000,5100] 

[5500,5600] 

[5600,5700] 

[5700,5800] 

[5800,5900] 

[5900,6000] 

[6000,6100] 

number of peaks 

20 

23 

24 

25 

25 

24 

23 

interval {N = 100000) 

[50000,50100] 

[50000,50200] 

[52000,52100] 

[54000,54100] 

[56000,56100] 

[58000,58100] 

[60000,60100] 

number of peaks 

20 

20 

21 

22 

24 

25 

23 


first the region close to a; = 0, Fig. [3[a) shows 15 peaks 
in the region [0, 300] for N = 1000. The corresponding 
region for N = 10000, which is [0,3000], contains 148 
peaks, confirming the conclusion drawn above that the 
frequencies are the same in corresponding areas for dif¬ 
fering N, leading to the overall linearity in N of the peak 
number (Fig. [8]). However, by counting the first 15 peaks 
in the probability distribution for N = 10000, shown in 
Fig. [7Kb), they fill the interval [0,960], indicating a \fN 
scaling of the maximum wavelengths around the center 
of the distribution. We clarify that this is not in con¬ 
tradiction with Table U where the representative low-x 
intervals are taken at different finite values of x. 

This type of scaling may also be observed in other re¬ 
gions of the QW probability distribution. Focusing next 
on the special structures noted above. Figs. [7Kc) and 
Figs.[7Kd) show three of these, which we find quite repro- 
ducibly in the region around x = N/2. More accurately, 
these macroscopic dips of the distribution envelope ap¬ 
pear around x = 0.36A^, x = 0.45A^, and x = 0.58A^. The 
interval around 0.58A^, which is also the region with max¬ 
imum oscillation frequency, shows a particularly remark¬ 
able beating structure [Figs.[7Ke) and[7Kf)], with multiple 
points where the upper and lower envelopes meet. Fig¬ 
ures [ZKg) and [ZKh) show the maximal frequency regime 
and parts of the beating envelopes in the fullest detail, 
and Table El contains the corresponding information for 
the beating interval for QWs of = 1000, 10000, and 
100000. Again it is clear that a factor-10 increase in N 
causes only a factor i/lO magnification of the width of the 
beating structure and, given the fixed maximal frequency 
in this interval, of the number of peaks it contains. 

We close our discussion of the scaling of peak widths 
by considering the probability oscillations around x = 
N / \/2. By counting the width of the region covered by 


TABLE 11. Position, width, and number of peaks in the dis¬ 
tinctive beating structure around x = 0.58A for walks of 
N = 1000, 10000, and 100000 steps. 


N 

region 

width 

number of peaks 

1000 

[546,604] 

58 

14 

10000 

[5682,5860] 

176 

44 

100000 

[57452,58016] 

566 

141 


the last 10 peaks up to and including the peak of max¬ 
imum probability, we find that this scales according to 
]Vi/3 (Fig. [^. Similarly, the full width at half maximum 
height (FWHM) of the leading peak in P{x) also scales 
with and hence this last peak retains its aspect 

ratio when x is normalized by N. Thus we demonstrate 
the presence of algebraic scaling in the probability oscil¬ 
lations over the full distribution. 


C. Peak Shape 

Next we consider the shape and scaling of the asym¬ 
metric envelope of peaks in the probability distribution to 
ascertain its functional form, P{x,N). The dependence 
on N is largely contained in Fig.|6|and we make these re¬ 
sults more systematic here. Concerning the dependence 
on x, we take the results of Sec. Ill A as a demonstration 
that the “crossover” region just beyond the maximum 
peak becomes a set of vanishing measure at large N-, to 
fit the envelope function, meaning the set of points ex¬ 
tracted from the full data set that fall close to the upper 
edge of Pix), for the region x < Xmax, we assume that 
it diverges at cc = ± Nfy/2 for large values of N. Before 
considering the QW, we recall that the functional form of 
the probability distribution for a classical random walk 
(Sec. IIA) becomes a Gaussian at large N, with the form 

P{x) = Po + Ae~ * 2 J , (9) 

where a = y/N, A = Ijy/2'KN, and Pq = 0 = 6 for a nor¬ 
malized and centered distribution. This is an exponential 
function whose characteristic width scales with '/N. 

By contrast, for the QW with large N we find an ex¬ 
cellent fit to an algebraic form, 

Pe{x) = Pq + a{b - xy (10) 

for each half of the distribution, with 5 = ± N/ y/2. To 
quantify the extent of the validity of such a fit, we exam¬ 
ine the data on logarithmic axes [inset. Fig. fTOT al]. find¬ 
ing that a single power-law provides a robust description 
of the envelope for the entire region 0.4A^ < x < a^max- 
For the three free parameters, our results as A^ —>■ 00 
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FIG. 7. Illustration of probability distributions in different 
parts of the interval for ID QWs of different step numbers 
N. First 15 oscillation peaks for (a) N = 1000 and (b) N = 
10000. Central region of the distribution for (c) N = 4000 
and (d) N = 10000, showing apparent beating. Detail of the 
envelopes in the region of the beating structure for (e) N = 
4000 and (f) N = 10000. Detail within the envelope of the 
beat structure for (g) N = 4000 and (h) N = 10000, showing 
the most rapid oscillations in the distribution. Probability 
maxima for (i) N — 100000 and (j) A'" = 1000000. 


(in practice, up to = 1000000) indicate that the con¬ 
stant Pq = —B/N with B = 1.884. Fits performed using 
only two remaining free parameters, and illustrated in 
the main panel of Fig. [TODa) for the case N = 1000000, 
allow us to deduce for the large-limit that the prefactor 
approaches a = A/\fN with A = 1.5 and the exponent 
approaches c = 0.5. We conclude that to an excellent ap¬ 


proximation the envelope function follows a square-root 
form in x measured away from x = ± N/^/2^ and this 
determines the algebraic form of the peak width. The 
behavior of the prefactor A ensures that Pe(x) oc 1/A^ 
across all of this range, consistent with the results of 
Fig.|6]but excluding the final peak. We remind the reader 
that the envelope function Pe(x) is not a quantity obey¬ 
ing a normalization law as the true distribution P{x) 
does. 

Below X = 0.4Af, the shape of the envelope begins 
to deviate from the universal square-root form [inset, 
Fig. linKa)]. To ascertain its shape close to the center 
of the distribution, we instead apply a fit of the form 

Pe(x)=Tb + aV' (11) 

for each half of the distribution. Again the data on 
logarithmic axes [inset. Fig. [rUT bl] show an excellent 
fit to a single set of parameters over a broad region, 
0 < X < 0.2N, where Pq = B'/N with B' = 0.615, 
a' = A' IN'^ with A' = 1, and c' = 2. The results are 
shown in the main panel of Fig. fTOT b'). Thus the P oc 1/A^ 
form is maintained, the universal behavior of the enve¬ 
lope around its center is a quadratic dependence on x, 
and there is a relatively significant constant contribution 
that reflects directly the incomplete nature of destructive 
interference in the central region of the QW. We regard 
the intermediate regime Q.2N < x < 0.4N as a crossover 
zone between the two limiting forms [Eqs. (fTUl) and dill) ] 
and do not consider it further. 

We close our discussion of the probability distributon 
P(x, N) by considering the region x > N/^/2. Flere there 
are no longer any oscillations [Figs. [71[i) and [71[j)] and 
P{x^ N) is the “envelope.” In Sec. IIIA [Fig. [6] we showed 
that the dependence on N crosses very rapidly to an ex¬ 
ponential decay around x = N/\f2. A complete fit of the 
data in this regime reveals the form 

P{x) = Pq + ae wu.B ^ (12) 



FIG. 8. Total number of peaks in the probability distribution 
for QWs from N = 100 to 100000. 






















































































FIG. 9. Width in x of different numbers of peaks in the prob¬ 
ability distribution for QWs over the full range of N values 
studied. The width of the first 10 peaks from a: = 0 (black) 
scales with y/N ; the width of the last 10 peaks up to the max¬ 
imum (red) and also the FWHM of the tallest peak (green) 
both scale with 

with Pq = 0, b = N/y/2, d = 3.2, and a = with 

A = 0.8, i.e. an exponential and pseudo-Gaussian behav¬ 
ior but with unconventional alterations to the exponents 
in both X and N. The effectiveness of this fit is shown in 
Fig. note), which also highlights how rapidly the proba¬ 
bility falls away in a short distance beyond x = N/y/2. 


D. Fourier Transformation of the ID-QW 

In Sec. IIIB we explored the rich spatial frequency 
information contained in the oscillations of the proba¬ 
bility distribution P{x). The QW contains oscillations 
on all length scales from ultra-short wavelengths around 
X = 0.58N to long-wavelength oscillations scaling as '/N 
around a; = 0, with similar algebraic scaling around 
X = Xmax- These differing spatial frequencies can even 
combine to create highly reproducible beating structures 
in certain regions. All of this information should be re¬ 
flected in the Fourier transform of P{x), which we dis¬ 
cussed from an analytical point of view in Sec. IIC. 

Here we perform a discrete Fourier transformation on 
the data sets for ID QWs of all lengths N, finding results 
for the Fourier components, F{k), of the type shown al¬ 
ready in Fig. |3Kb) for N = 100 and N = 1000. F{k) 
possesses a primary peak with amplitude F{0) = 1 at 
fc = 0, flanked by two secondary peaks with negative 
components at fc ~ and then shows an oscillatory 
form between positive and negative values of the Fourier 
components. The oscillations are again contained within 
a decaying envelope function, which we find to be iden¬ 
tical at positive and negative values, and the fc-space pe¬ 
riodicity of the oscillation is remarkably constant across 
the range — tt < A: < tt. This result is a clear reflection of 
the fact that spatial information is present in the prob- 
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FIG. 10. Fit to the envelope of the probability distribution, 
Pe{x), for a QW of A = 1000000 steps, (a) Fit for the interval 
400000 < X < XmBx using the algebraic form of Eq. m with 
B = 1.884, A = 1.487, b = 707144, and c = 0.4918 (see 
text), (b) Fit for the interval 0 < a; < 200000 using the 
algebraic form of Eq. (Hill with B' = 0.6152, A' = 1.0017, and 
c' = 2.074 (see text), (c) Fit to the probability distribution 
P{x) for the interval x > Njypi using the exponential form 
of Eq. (fT21l with Pq = 0, b = N/ y/2, d = 3.2, and A = 0.8 (see 
text). 


ability distribution on all length scales, and the mixture 
of positive and negative components across the range of 
k is manifest in complex mixing phenomena such as the 
beating structure. However, beyond the large fc = 0 com¬ 
ponent there are no special spatial frequencies appearing 
in the distribution. We comment here that the constant 
component F'(O) = 1 is simply the sum of all data in real 
space, and therefore is the result expected for a normal¬ 
ized distribution. 

We begin our quantitative analysis of the Fourier trans¬ 
form F{k) by counting its peaks. Figure ITlT a) confirms 
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FIG. 11. (color online) (a) Number of peaks in Fourier space 
for QWs of lengths N = 100 to 100000. (b) Length of the 
envelope function of the final beat below fc = tt, in units of 
Nk/2n (black), and the number of peaks in this region (red), 
shown as a function of y/N. 


that the total number of peaks in the Fourier spectrum 
scales linearly with TV, as in real space and again with a 
constant of proportionality of order 8% (more precisely, 
1/12) for each half of the transformed distribution. As 
in real space, we may also count the peaks in particu¬ 
lar parts of the distribution to investigate their scaling 
form. The most striking feature of the Fourier trans¬ 
formed data for large values of is a reappearance of 
beating phenomena between the upper and lower enve¬ 
lope functions. In complete consistency with the results 
of Sec. IIIB, where the beating structures were observed 
in the region of the distribution with the highest spatial 
frequencies, the Fourier-space beats are clearest close to 
k = TT. In Fig. [T^ we illustrate this property with the 
fc-axis rescaled to Nk/2TT to better reflect the number 
of Fourier components in the data set. As in Sec. IIIB, 
we may characterize the structure of the beat pattern by 
considering the length of the final beat and the number 
of peaks it contains, which are tabulated in Table Hill and 
illustrated in Fig. ITlT b). From the latter it is clear once 
again that the beat structures scale according to y/N. 

Proceeding as in Sec. IIIC, we consider the possibil¬ 
ity of a universal fit to the envelope function Ff.{k,N), 
whose fc-dependence is valid for all values of N. The N- 
dependence of Fe{k, N) is shown in Fig. ITSl a) for values 
of k from across the full range, and the constant slope 
gives the clear result Fe{k^N) oc l/y/lv. Mindful of the 
fact that the real-space envelope changes form between 
the limits of small and large x, we consider the functional 


TABLE III. Characterization of oscillation frequencies in 
Fourier space. The length of the last beat is quoted in units 
where the interval of the Fourier components is rescaled to 
{-Nk/2n,Nk/2Tv]. 


N 

total peaks 

peaks in last beat 

length of last beat 

1000 

167 

9 

45 

4000 

667 

18 

87 

10000 

1667 

29 

140 

40000 

6667 

58 

276 



FIG. 12. Fourier components of the probability distribution 
close the maximum frequency, shown rescaled to Nk/2'K. In 
this regime the envelope function shows clear beating behav¬ 
ior. 


forms 


and 






(13) 


(14) 


for small and large k. As shown in Figs. fTSl b) and fTSl c). 
Eq. (1131) with exponent c = 1/2 provides an excellent fit 
at small k, out to /c ~ 0.27r, and Eq. m with c' = 1 an 
excellent fit for all k values in the upper half of the range. 
Thus we find the fc-dependence of the envelope of Fourier 
components to be algebraic over the whole range, with 
an inverse square-root decay away from fc = 0 crossing 
over to a linear decrease as fc approaches tt. 

In summary, the Fourier transform of the QW prob¬ 
ability distribution contains all of the same information 
in a complementary form. It is bounded by upper and 
lower envelope functions with the same algebraic decay. 
It demonstrates that spatial frequencies are present on 
all scales from the inverse step length to the inverse 
system size, with no special internal period(s) but with 
distinctive beating structures on a length scale of y/N. 
While the QW does not satisfy the strict definition of self¬ 
similarity (no fractal structures appear), many of its fea¬ 
tures are similar and scale-invariant across the full range 
of N values. Thus the simple quantum evolution algo¬ 
rithm of Sec. IIB contains a very rich variety of spatial 
information. 
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FIG. 13. Fit to the envelope of the Fourier transform of 
the probability distribution, Fe{k), for QWs of = 1000 
to 100000 steps, (a) A^-dependence of F^{k) for four different 
values of k, showing a clear algebraic decay with exponent 1/2. 
(b) Fit using the algebraic form of Eq. (1 1311 . demonstrating its 
relevance for the interval 0 < k < 0.27r and returning the ex¬ 
ponent c = 1/2. (c) Fit using the algebraic form of Eq. (I14II . 
demonstrating its relevance for the interval OAtt < k < tt and 
returning the exponent c' = 1. 


IV. TWO-DIMENSIONAL QUANTUM WALKS 
A. Classical Random Walk in 2D 

As noted in Sec. I, it is the advent of experiments capa¬ 
ble of realizing a 2D QW [2^ that has caused the resur¬ 
gence of interest, both experimental and theoretical, in 
the field. Before discussing the range of properties exhib¬ 
ited by QWs in 2D, it is helpful to review the classical 
case. In principle there are two ways to generalize the 
unbiased ID classical random walk to 2D. 


1) Adopting a square grid, the walker is equipped only 
with a two-face coin and therefore flips it once to step to 
(x ± 1, y), then a second time to arrive at (x ± 1, j/ ± 1). 
Now the conventional V-step ID binomial distribution 
for a: or 2 / is recovered by summing over the probabilities 
in the orthogonal direction and it is easy to show for any 
N that P{x,y) = P{x)P(y). In the large-A limit, the 
distribution approaches 

P{x,y) = (15) 

and the mean distance of the walker from the origin after 
N steps is (r) = \/l^. 

2) Remaining on a square grid, the walker at site (x, y) 
has equal probabilities of 1/4 (equivalent to a four-face 
coin) to move to any of the points (x± 1, j/i 1). It results 
in an identical distribution to case (1). 

A further generalization is to consider a continuous 
space in which the walker takes unit-length steps at any 
random angle 0 < (j) < 2 tt. The walk approaches perfect 
circular symmetry (r^ = + y^) at large values of N. 

B. Two Types of 2D QW 

With a view to experimental realization, 2D QW also 
has two-face coin scheme and four-face coin scheme. 

A four-face coin can be constructed by product 
the Hadamard operators, i.e. H 4 = Uh (S' Uh{b.s well as 
corresponding initial state). It evolves to a distribution 
shown in Fig. [14] for a walk of 100 steps. We note imme¬ 
diately that this 2D QW is a directly expand of ID QW, 
in which P{x,y) = P{x)P{y). One will understand this 
equivalence by invoking the matrix identity 

{AiBi) 0 {A 2 B 2 ) = {Ai 0 A2 ){Bi 0 B 2 ), (16) 

Ai , A 2 represent the Hadamard operator Uh and Bi , B 2 
for the initial state | t) + *1 i)- 

The two-face-cion scheme lead to a completely differ¬ 
ent 2D distribution, shown in Fig. [15] for a walk of 100 
cycles. Specifically, an evolution protocol using the same 
two-face coin twice in each complete cycle to generate 
successive steps in the x and y directions. The first flip 
of the coin selects the direction of ±x and the second zty. 
A full cycle can be represented as tj = SyY(SxY) and 
the quantum state after N evolution cycles as I’ipw) = 
(17)^|'0o)-This probability distribution has several prop¬ 
erties in common with the ID QW, including strong prob¬ 
ability peaks far from the center, strong destructive in¬ 
terference everywhere near (x, y) = (0,0), and oscillatory 
behavior around the locus of maxima. However, the qual¬ 
itative features of this QW are strikingly different from 
those of the four-face-coin walk, most notably in that 
the peaks are rotated by 45°, they occur right at the 
edge of the system, and here the locus of maxima does 
show some tendencies towards a circular shape, even if 
the peaks upon it remain strongly anisotropic and four¬ 
fold symmetric. The location of the probability maxima 
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is really very anomalous, in that the minimum of the 
destructive interference is truly obtained when a walker 
makes all N of its steps in the same direction for one lat¬ 
tice orientation, but precisely N/2 in each direction for 
the other orientation. 



FIG. 14. Probability distribution of 2D QW effected using 
a four-face coin selecting steps in the ±(x-|-j/) and ±{x-y) 
directions, for N = 100 steps. 



FIG. 15. Probability distribution of 2D QW effected using a 
two-face coin twice in each cycle to select steps in the ±1 and 
±y directions, for N = 100 cycles. 


Fi g, fibl can also be realized by a four-face coin(Grove 
walk|ll|b In one of the most notable recent studies [H, 
[T^ . Di Franco and coauthors proved the equivalence of 
the spatial probability distributions between the 2D QW 
using a single two-face coin, which these authors termed 
the “alternate quantum walk” (AQW) and the Grover 


walk. The AQW is actually a 2D QW in which the two 
direction are maximal entangled. 

The Grover coin can be expressed in the form 


G4 


1 

2 


/-I 

I 
I 

\ 1 


1 

-1 

I 

I 


1 

I 

-1 

I 



(17) 


Noted that the distribution and spreading rates are ef¬ 
fected by the initial state fill , fl^ . EqUH] is the one pro¬ 
viding Grove walk a maximum spread rate 

IV'™) = id t) - I -t) - I f-) + I i))|0x,0j;). (18) 

We conclude here the difference between the two type 
of 2D distribution, is not distinguished by two-face coin 
or four-face coin, but whether the two orthogonal direc¬ 
tions are entangled. Fig.[T3]can also be realized using two 
uncorrelated two-face coins. Some authors have studied 
decoherence [1^ and localization [s^ in 2D QWs using 
the concept of two coins, and others have quantified the 
effects of entanglement by a partial or complete swapping 
of the two coins after each step of the walk 

The method of generating a QW by using the same 
two-face coin twice, to which we refer henceforth as the 
AQW P,[l3, is an important one for a number of rea¬ 
sons. First and foremost is that two-state systems are 
much easier to find, or to produce, in any physical real¬ 
ization of a quantum walker, and hence are much more 
relevant for experimental implementation. Secondly, this 
maximally entangled protocol contains further unconven¬ 
tional phenomena, which we discuss in Sec. V. Further, 
coin entanglement provides a valuable and completely 
general means of constructing and perhaps also of real¬ 
izing high-dimensional QWs using only two-face coins. 

We conclude this subsection by contrasting the two 
QWs generated by unentangled and by maximally en¬ 
tangled coins. The unentangled system appears to show 
a perfectly square, factorized probability distribution 
(a result we demonstrate in Sec. V) characteristic of 
the ID QW, with maxima at positions converging to 
(± A^/-\/2, ± Af/-\/2). The entangled case shows a circu¬ 
lar “probability front” and significant complexity in the 
interference pattern within it. In the sense that the un¬ 
entangled situation can be discussed as accelerated diffu¬ 
sion from the classical random walk to the quantum walk 
as a result of destructive interference, so the maximally 
entangled walk appears to show a still further accelerated 
diffusion. Indeed, it achieves a situation where informa¬ 
tion propagates to the very edge of the system with high 
probability, which is a quite remarkable consequence of 
near-perfect destructive interference among trajectories 
in the center region. A possible interpretation of this 
result may be found in preservation of the information 
content of the coin, because the degree of two-coin en¬ 
tanglement is “complete” in the sense that the walk can 
be generated using only one coin. 
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V. NUMERICAL STUDIES OF 2D QWS 

For an exact characterization of the 2D QWs intro¬ 
duced in Sec. IVB, we turn now to a numerical investiga¬ 
tion of their probability distributions. We begin by con¬ 
trasting in Table HVl the analytical results for the proba¬ 
bilities of the two walks for = 6 steps. 

This exercise makes clear that the two QWs are radi¬ 
cally different from the outset. For the unentangled case, 
it is easy to see the result one may already suspect from 
Fig. [T31 that the probability distribution of the unentan¬ 
gled 2D QW, P{x, y) = P{x)P{y), is a direct product of 
two ID QWs in the x and y directions. We will shortly 
demonstrate it numerically for large values of N. We 
remind the reader that this result is not exactly intu¬ 
itive, given that the walk steps are each made in the 
x±y direction of the lattice and as such would appear 
to entangle the two directions completely; however, this 
result also emerges from the 2D classical random walk 
(Sec. IVA). For the entangled case, the probability dis¬ 
tribution has no direct relation to the ID QW and indeed 
already shows the concentration of probability along the 
X and y directions rather than along the diagonals. 


A. Cross-Sections of the Probability Distribution 

For our numerical calculation of the characteristics of 
the different 2D QWs, we have computed the probability 
distributions for both walks up to iV = 1000 steps, whose 
illustration requires a 1000x1000 grid. To show the re¬ 
sults in a manner more quantitative than is possible in 
Figs.fMlandlTSl we consider slices through the 2D prob¬ 
ability data taken along the x axis, along the diagonal 
x+y, and along the “edge” of the data set, as shown in 
Fig. [m 

The results of this process are illustrated in Fig. [T7] 
for walks of iV = 100 steps. For the non-entangled case. 


TABLE IV. Probability distributions of the 2D QW for A = 6 
steps, shown as 4096P(a;, y). On the left is the four-face- 
coin QW (equivalent to two unentangled coins) and on the 
right the AQW (equivalent to a single coin, or two maximally 
entangled coins). For the interpretation of these numbers we 
note that in the ID QW with A = 6 one obtains 64P(a;) = 
[1,18,9,8,9,18,1]. 
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FIG. 16. Schematic illustration of ID slices taken through the 
2D probability distribution for the two 2D QWs considered, 
whose data are shown in Fig. 1171 Slice A denotes the x- 
direction, B the diagonal, and C is the edge of data. In panel 
(a), where the blue square denotes the region of maximum 
probability, slice C is taken at a; = 0.7A; in panel (b), where 
the maximum probability is found along the blue circle, slice 
C is taken at the true edge {x = A). 


panel A1 proves the numerical identity with the ID QW, 
P{x, 0) = P{x)P{0), which can be compared with Fig. [5] 
using the result for P(0) with N = 100. This being the 
case, panel Cl is very easy to interpret and for any par¬ 
allel cut would have the same functional form with any 
prefactor from P{x). Panel B1 can be expected to satisfy 
P{x, x) = P^{x), and therefore to have an envelope func¬ 
tion of the form {N/\/2—x)~^ across the outer half of the 
distribution, a result we demonstrate below. By contrast, 
rather little is known about the distribution P{x, y) for 
the AQW and its understanding will require applying the 
techniques of Sec. Ill, for small data sets, to the panels 
A2, B2, and C2 of Fig. [T6l Qualitatively, in A2, the hor¬ 
izontal slice through the probability maxima, we observe 
an apparent envelope function with no oscillations and a 
divergence towards x = ±A with an unknown power. In 
B2 we observe a complex oscillatory pattern with signifi¬ 
cant probability extending beyond its peak value. In C2 
we observe a remarkably classical-looking probability dis¬ 
tribution at the edge of the system, where we remind the 
reader that the peak is the absolute maximum anywhere 
in the walk. 

We wish to characterize the AQW by the power-law 
form of its probability envelope function around the peak 
values of the distribution. For this we analyze the 2D 
probability slices following Eq. OT, but because our 
largest system size in 2D is A = 1000, the numerical re¬ 
sults do not give particularly reliable curve fits. Although 
the divergence of P(x) a.i x = Nl\/2 (at a; = A for slice 
A2) is not achieved with any accuracy for these values 
of A, we enforce this value in all cases other than slice 
B2 to reduce the arbitrariness in the fitting parameters. 
We show in Fig. [TH] the probability on logarithmic axes 
as a function of (6 — x)/N for fixed values b — 0.707A in 
cases A1 and Bl, 6 = A in case A2, and the fitted value 
b = 0.8A in case B2. As expected, the probabilities for 
panels Al and B1 have gradients close to —1/2 and —1 
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FIG. 17. Probability distributions on ID slices through the 
2D QW data (see Fig. fThll for walks of = 100. Panels Al, 
Bl, and Cl are for the unentangled QW and A2, B2, and C2 
for the AQW. 


FIG. 18. (color online) Scaling properties of the envelope 
function of the 2D probability distribution for N = 500 
(black) and N = 1000 (red) on the data slices shown in 
Fig. 1171 panels Al, Bl, A2, and B2. 

although we point out that for x ^ N/ v^, where the ID 
QW shows its strongest peak, this distribution has its 
deepest valley. 


[c = 1/2 and c = 1 in Eq. (|T0|) ]. respectively, over the bulk 
of the range, and the accuracy to which the data fall on a 
straight line benchmarks the method for N = 1000. For 
panel A2, the data are remarkably clean and give a clear 
gradient parameter c = 1, indicating that the envelope 
function of the peak in the maximally entangled walk 
satisifes the form P{x,0) = Pq + a/{l — x/N) to high 
accuracy. For panel B2, a slice that does not include the 
main peaks but only the circular edge of local probability 
maxima, we find a result close to c = 1/2, although here 
the envelope is poorly dehned and the errors significant. 

While our estimate of the functional form of the x- 
dependence of the AQW is somewhat approximate, the 
^dependence of the probability distribution is beyond 
doubt. At every corresponding position, the probability 
falls by a factor of 4 for every doubling of N. This result 
is illustrated by the black (A^ = 500) and red {N = 1000) 
lines in Fig. [TS] Thus we propose that the appropriate 
fit to the probability data on these four slices is given by 
the algebraic form 


y = ^ + (19) 

with c+d = 2 and the values of c as deduced from Fig.jTH] 
We may conclude that the entangled 2D QW has alge¬ 
braic scaling properties similar to the ID-case. While its 
probability distribution shows no oscillations in the x and 
y directions, in the x±y directions it oscillates in a man¬ 
ner not dissimilar to the ID QW (panel B2 in Fig. fTTl) . In 
Fig. [THl we show that the number of peaks in this slice is 
again proportional to the step number N, and with a con¬ 
stant of proportionality again of order 8.5% (Sec. IIIB), 


B. Edge of the Probability Distribution 


Next we consider the probability distribution at the 
edge of the walk for the AQW, shown in panel C2 of 
Fig. [T71 As remarked above, this QW has the highly 
anomalous feature that its maximum probability occurs 
when the walker takes the maximum number of steps in 
the same direction along one of the two axes, but returns 
to the center of the other axis. Further, a visual inspec¬ 
tion of the ID slice through this maximum suggests that 
the distribution on this edge may be a Gaussian, or a 
related function similar to the ID classical random walk. 

To investigate whether a quantum walk can lead to 
a classical results, we employ the decomposition of the 
unitary evolution matrix into the matrices P and Q in¬ 
troduced in Sec. IIIBI For the 2D QW we require ma¬ 
trices Px for steps to the left, Qx for right, Py for up, 
and Qy for down. Exploiting the AQW equivalence of 
the maximally entangled two-coin QW, we use the same 
coin alternately for the x and y directions. Because our 
aim is to understand the probability distribution on the 
right-hand edge of the system (equivalent to Fig. fTHT b). 
slice C), every ^-direction operation for the walker must 
be a Qx matrix and not Px ■ A complete cycle of the walk 
may then be either R — PyQx or S = QyQx, where 


R = PyQx 
S — QyQx 


1 

2 

1 

2 




1 fo 0 


2 Vl -1 


( 20 ) 
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FIG. 19. Number of peaks in the probability distribution of 
the AQW (maximally entangled two-coin QW) in the direc¬ 
tion x+y, shown in panel B2 of Fig. 1171 as a function of step 
number N. 
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FIG. 20. Comparison between the AQW edge probability dis¬ 
tribution and a binomial (classical random walk) distribution 
for values up to A = 1000 steps, where the binomial is very 
accurately Gaussian, (a) Standard deviation a. (b) Prefactor 
A. 


RSRR, and SRRR, yielding 


^ RRSR + RSRR = ^ ^ 


1610 or 
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and thus 


p{N,2) _ p(JV,2) p{N,2) _ foal 

^ - 256' ^ ’ 


The only path arriving at position {N, 4) is RRRR, lead¬ 
ing to 


i.e. only the left-most operator in the string of steps is 
important. 

With this result it is possible to calculate the entire 
edge distribution analytically for any value of N. We il¬ 
lustrate the process for a walk of A = 4 steps. The paths 
arriving at position (A, 0) are RRSS, RSSR, RSRS, 
SSRR, SRSR, and SRRS, which are divided between 
the coin states | f) and | 1) according to 

= RRSS + RSSR + RSRS = 

= SSRR + SRSR + SRRS = 


1 /3 -3 


16 VO 0 

1 /3 -3 


16 VO 0 


whence 
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( 22 ) 


Paths arriving at position (A, 2) are RRRS, RRSR, 


= RRRR = 
p{NA) _ piNA) _ 

^ ~ rt) - 


16 Vo 

1 — i 


16^2 



(24) 


This analytical solution demonstrates that the quantum 
mechanical interference leading to the probability on the 
edge of the system is constructive everywhere, with none 
of the paths interfering destructively. Regions at the cen¬ 
ter of the edge simply have the most paths, and this is 
the origin of what we can call the “semi-classical” result 
that the probability is maximal at the center of the edge, 
leading to a distribution similar in appearance to the bi¬ 
nomial distribution of the classical random walk. In fact 
the degree of destructive interference everywhere else in 
the 2D AQW is such that these maxima on the edges are 
the global maxima, meaning that the walker is effectively 
pushed to the maximal number of steps in order that it 
does not “destroy itself” by interference in the maximally 
entangled two-coin walk. 

A quantitative examination of the coefficients of the 
AQW edge reveals that, despite being peaked at the cen¬ 
ter, they are not the same as the ID classical random 
walk, as shown in Table IS The probability distribution 
at position x on the edge of an A-step AQW can in fact 





















15 


be expressed exactly as 


Pn{N,x) 


~ 1 ) + Pn-i{^ + 1 )] 

1 

4IV 


C 


N-l 


+ 


c 


N-l 


(25) 


where is a binomial coefficient (Sec. IIA). This re¬ 
markable result, which we term a “pseudobinomial” dis¬ 
tribution, encodes the emergence of the physics of clas¬ 
sical probabilities in the highly entangled QW. The bi¬ 
nomial coefficients appearing in the penultimate step be¬ 
fore measurement (Table |V]) arise as a consequence of 
the complete constructive interference of paths, as illus¬ 
trated in Eqs. (|22l) to (|24ll and discussed in the preceding 
paragraph. 

To demonstrate the form of the AQW edge probability 
at large N, we fit both binomial [Eq. ([1])] and pseudobino¬ 
mial distributions to the Gaussian, as specified in Eq. 
and again with Pp = 0 = 6, to compare their forms 
and to extract their standard deviations a. As shown 
in Fig. EUT al. while the binomial distribution gives the 
result cr = \/]V, the st andar d deviation of the AQW edge 
probability is ct = ^yN/2, i.e. the distribution is nar¬ 
rower and the diffusion, or spreading, rate of the walk 
slower (more localized) by a factor of l/-\/2. Concerning 
the normalization prefactor A [Fig. EUT bl]. the binomial 
approaches the Gaussian result A = but the 

AQW edge distribution does not, following instead a de¬ 
pendence A oc 1/N. While the probability of the classical 
random walk of course sums to unity, determining A for a 
true Gaussian, the probability distribution at the edge of 
the AQW is not normalized due to the weight in the inte¬ 
rior of the walk, and hence the edge probability is found 


TABLE V. Probabilility tables for the first 7 steps of the bino¬ 
mial distribution (top), omitting a prefactor of 1/2^ at each 
level N of the table, and for the pseudobinomial distribution 
(bottom) achieved on the edges of the 2D AQW, omitting a 
prefactor of 1/4^. 
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FIG. 21. 2D Fourier transform F{kx,ky) of the probability 
distributions P{x,y) (Fig. ??) for (a) the four-face-coin or 
unentangled 2D QW and (b) the 2D AQW or maximally en¬ 
tangled two-coin QW, both for N — 100. 
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to decay more rapidly than a true Gaussian. The emer¬ 
gence of such a pseudobinomial distribution in a QW is 
yet another example of the rich physics contained in a 
deceptively simple quantum evolution algorithm. 


C. Fourier Transformation of 2D QWs 

For further insight into the structure of the AQW prob¬ 
ability distribution, we compute the Fourier transform 
F{kx,ky) of P{x,y) for both the 2D QWs we consider. 
Figure I^TT a) shows the 2D Fourier transform of the un¬ 
entangled (four-face-coin) QW, whose product structure 
is again clearly visible. It was shown in Sec. Ill that the 
ID QW has spatial frequency information at all scales up 
to the inverse step size {k = tt) and this is clear in the 
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FIG. 22. Cross-sections of the 2D Fourier transform F{kx, ky) 
of the probability distribntions [Fig. [T7] for the two 2D QWs 
of Sec. IV, calculated with N = 100 and shown for the slices 
stndied in Sec. VA fFig. ITOl). 


finite components oi F{kx,ky) up to the edge of Fourier 
space. The AQW, shown in Fig. [2T][b), is again quite 
different, showing both an intrinsically 2D character and 
an apparent concentration of weight closer to the center 
of Fourier space. 

To analyze these distributions in a quantitative man¬ 
ner, we present the data in the form of ID slices. Figure 
[22] is completely analogous to Fig. |T3 with the unen¬ 
tangled QW in the left panels, the entangled one on the 
right, A denoting a horizontal slice through the center of 
the Fourier distribution, B a diagonal slice, and C the 
edge. We illustrate the qualitative features of the data 
using QWs of = 100 steps, but for our numerical anal¬ 
ysis of the properties of the Fourier transforms we use 
values of N up to 1000. 

In panel AI of Fig. [211 which shows F{kx,0) for 
the unentangled QW, we observe the factorized form 
F{kx)F{0) = F{kx) expected from Sec. VA and shown in 
Fig.EIb). Panel BI shows F{kx, kx) = F'^{kx), with only 
positive components and a correspondingly steeper decay 
of the envelope. Panel Cl, showing F{TT,ky), is identi¬ 
cal to Al up to a multiplicative prefactor, which is small 
and happens to be negative at kx = tt. Turning to the 
AQW, where the probability distribution cannot be fac¬ 
torized, the slice F{kx,0) in panel A2 is not dissimilar to 
panel Bl, in that all of the Fourier components are posi¬ 
tive and they decay significantly more rapidly than those 
of the ID QW. However, their oscillatory form shows a 
very precise odd/even modulation, which is not evident 


in the ID QW. The diagonal slice F{kx, kx) is shown in 
panel B2 and confirms both the odd/even modulation 
and the concentration of Fourier amplitude near the cen¬ 
ter of the system. Finally, panel C2 has no readily dis¬ 
cernible structure, reflecting the fact that the AQW edge 
is non-oscillatory and thus dramatically different from all 
the other QW distribution slices (Sec. VB). 

Following Sec. HID, we investigate two properties of 
the Fourier transform slices, namely their envelope and 
their oscillations. To characterize the decay of the enve¬ 
lope function, we follow the procedure shown in Fig. 1131 
and use logarithmic axes. As shown in Fig. [^5f ab and as 
expected from Sec. HID, the envelope functions for pan¬ 
els Al and Bl yield algebraic decay exponents close to 
c = 1/2 and c = 1, respectively, for the Fourier compo¬ 
nents around fc = 0; the accuracy with which the data ad¬ 
here to a straight line again benchmarks the accuracy one 
may expect from data up to V = 1000. The data from 
slice A2 form by far the best-quality set in Fig. [^5f al and 
the algebraic decay exponent is unequivocally c = 1/2, 
meaning that the spatial periodicity content of the AQW 
is qualitatively similar to that of the unentangled QW 
at low frequencies. We take this opportunity to remind 
the reader that there is a priori no direct connection be¬ 
tween the real- and reciprocal-space exponents of the en¬ 
velope functions (respectively around the peak of P{x) 
and around k = 0) and the two QWs present an exam¬ 
ple pair with different exponents in real space, describ¬ 
ing the peak shape, but the same exponent in reciprocal 
space, describing the frequency content. The data from 
slice B2 exhibit the lowest-quality envelope in the figure, 
but still show a strong qualititive similarity to slice Bl, 
with c = 1. Again this highlights the complexity of the 
spatial frequency content of the AQW in its different di¬ 
rections and indicates that the entanglement inherent in 
the AQW also entangles the spatial information of the 
different lattice directions. 

For completeness, we show in Fig. [^5f bl the Fourier 
components of the four slices at high fc, where once again 
slices Al and Bl, with respective gradients c' = 1 and 
F = 2 benchmark the accuracy of the approach. Again 
we observe that slices A2 and B2 for the AQW are qual¬ 
itatively similar, with the caveat that the data for slice 
A2 show evidence of a distinctively different intermedi¬ 
ate regime. We stress that the difference in amplitude 
of the Fourier components between the unentangled QW 
and the AQW, which reaches two orders of magnitude 
over the large-fc half of the range for slices Al and A2, 
is the most meaningful measure of strong quantitative 
differences between the two QWs. 

Concerning the oscillations visible in the Fourier trans¬ 
forms, in the unentangled case it is known from Sec. Ill 
that the 2D QW contains spatial frequency information 
on all length scales and that beating phenomena become 
visible both in real space around the maximum frequency 
and as a result in Fourier space near k = ±7r. For 
the AQW, our data (Figs. and |23|) show that the 
high-frequency spatial components are very significantly 
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and deduce the eigenvalues 

= ± 1 , 

= -A[cos(fci + fe) + cos(A;i - fe) (26) 

± i/(cos(fci + k 2 ) + cos(A;i - k 2 )Y - 4]. 


Because Xk = the dispersion relations for the four 
eigenmodes of the AQW take the form 


1+ r, 1- 

ujI^ = 7r=F i(cosfci +cosA: 2 ), 


(27) 



N(1-k/7r)/2 


FIG. 23. Fourier transform of the 2D probability distribu¬ 
tions for slices Al, Bl, A2, and B2 shown in Fig. I17I (a) as 
a function of k, showing algebraic decay of the small-fc com¬ 
ponents, and (b) as a function of 1 — fc/rr, showing different 
power-law dependences at large k. 


weaker, meaning that lattice-scale oscillations are of little 
relevance, and this suggests that beating (which we are 
unable to find up to = 1000) is unlikely to be present. 
We believe that the extra strength of even harmonics in 
the AQW results from the even step number. 

We close our analysis of the Fourier-space information 
contained in a QW by summarizing a complementary 
perspective, named the “dispersion relation” approach 
in Ref. M- In Sec. lie we presented the exact solution 
of the ID QW by considering its nature in Fourier space, 
finding that the Hadamard operator ([3]) gave rise to a 
dispersion relation for a quantity oJk = sin~^[sinA:/\/2], 
specified in terms of k. Although ujk is largely a simplify¬ 
ing notation, as there is no concept of a separate “time” 
(step number) and space in a QW, it does result in a 
compact description in higher dimensions and it has the 
added advantage of reflecting different degrees of coin en¬ 
tanglement in a transparent way. When the AQW in 2D 
is transformed into Fourier space , one may express 
the operator as 



/ _g*fclgifc2 

gifeigifc2 

gifeigife2 

gifei gife2 

1 

gifei g-*fc2 

T 

1 

gi/cig-ifc2 

gifeig-ifc2 

2 

g-ifei gife2 

g-ifcigi/c2 

_g-*felg*fc2 

g-ifeigife2 


1 g-ifclg-ife2 

g—ifcig-ifc2 

g-ifcig-ifc2 

_e-ikie-i 


allowing a considerable simplification of the process 
[Eq. ([7])] of generating the probability distribution in this 
case. 


VI. SUMMARY 

Quantum walks produce probability distributions en¬ 
tirely different from the well-known classical “drunk¬ 
ard’s walk.” In most cases, the distribution is controlled 
completely by predominantly destructive interference be¬ 
tween the paths returning to the center, with the result 
that the positions of maximal probability are pushed out 
towards the edges of the walk. This results in the “lin¬ 
ear diffusion” property of a quantum walker, which has 
made the QW evolution algorithm so valuable in studies 
of quantum computing. 

However, as quantum computing approaches large- 
scale implementation, it is necessary to understand the 
nature of quantum algorithms at large scales. By this is 
meant an exact account of their properties, information 
content, limiting behavior, and scaling characteristics at 
long evolution times. Although analytical studies of dif¬ 
ferent quantum walks in one and two dimensions have 
revealed certain characteristics of their evolution and in¬ 
terference, in particular their dependence on the initial 
state and the entanglement of the quantum coins gener¬ 
ating them, to date there has been rather little consider¬ 
ation of the situation at large step number N. 

By calculating the probability distributions for one- 
and two-dimensional quantum walks up to iV = 1000000 
in the former case and N = 1000 in the latter, we have 
revealed a number of properties and scaling character¬ 
istics. In one dimension, we verify that the probability 
approaches peaks at x/N = ±l/v^ at large N, which 
is the transition region between destructive interference 
and vanishing probability for all steps to be made in the 
same direction. The normalized support converges to 
the region [— 1 /\/2,1 /\/2] and the envelope of the distri¬ 
bution peaks has a square-root decay, i.e. an algebraic 
form. Within this envelope, the probability shows sys¬ 
tematic oscillations on all length scales, with the number 
of probability peaks always the same fraction of N. The 
frequency of these oscillations increases from the inverse 
system size to the inverse step length as a function of 
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distance from the center of the walk. The different fre¬ 
quencies show complex beating phenomena in the regime 
where the oscillations are most rapid. These properties 
are revealed in a complementary fashion by taking the 
spatial Fourier transform of the distribution. All of these 
features are universal for walks of all N values, giving 
them very strong similarities, but not the property of 
self-similarity (there are no fractal structures in the sim¬ 
ple walks studied here). 

In two dimensions there is not one quantum walk, or 
even one algorithm, but a spectrum of protocols capa¬ 
ble of generating quantum evolution in a plane of phase 
space. We study two examples that in fact represent the 
limiting cases of the entanglement between two orthogo¬ 
nal directions: unentangled 2D QW and the maximally 
entangled alternate quantum walk(AQW). 

In our numerical studies of these two limits, the under¬ 
standing developed in one dimension gives a complete 
account of the unentangled quantum walk, whose two- 
dimensional probability distribution factorizes exactly 
into two one-dimensional functions. By contrast, the 
maximally entangled case exhibits strong correlations be¬ 
tween the two orthogonal directions, damping of the os¬ 
cillatory behavior, and the extraordinary feature that the 
maximum of the probability distribution is pushed all the 
way to the system edge by the dominance of destructive 


interference. We provide an analytical description of the 
edge distribution, showing that all paths arriving at the 
system edges interfere constructively and proving that 
its functional form is a type of pseudo-binomial, which 
is semi-classical in the sense of approaching a Gaussian 
dependence on the spatial coordinate at large N. 

As two-dimensional quantum walks become an exper¬ 
imental science, our analytical and numerical studies 
demonstrate that even the simplest algorithms for quan¬ 
tum evolution contain a rich variety of physical phenom¬ 
ena and potential for technological application. 
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